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For the spinor Bose-Einstein condensates both the total spin S and its Z-component Sz 
should be conserved. However, in existing theories, only the conservation of Sz has been taken 
into account. To remedy, this paper is the first attempt to take the conservation of both S and 
Sz into account. For this purpose, a total spin-state with the good quantum numbers S and Sz 
is introduced in the trial wave function, thereby a generalized Gross-Pitaevskii equation has been 
derived. With this new equation, the ground bands of the '^'^Na and *^Rb condensates have been 
studied, where the levels distinct in S split. It was found that the level density is extremely dense in 
the bottom of the ground band of ^^Na, i.e., in the vicinity of the ground state. On the contrary, 
for *^Rb, the levels are extremely dense in the top of the ground band, 
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1, Introduction 

In 1998 the spinor Bose-Einstein condensation of 
sodium atoms was experimentally realized in an optical 
trap by Stamper-Kurn et al 0. The advantage of the 
optical trap is that it liberates the spin degrees of free- 
dom of atoms (e.g., ^^Na, ■^^K, and ®^Rb atoms have 
nuclear spin 3/2 and electrons at s orbits, they behave 
as simple bosons with spin one at low temperature). As a 
consequence of the liberation, such a macroscopically co- 
herent system with internal degrees of freedom becomes 
very rich in physics. Several interesting phenomena have 
already been obtained, e.g., quantum entanglement of 
spins 0, 13 , spinor four- wave mixing 0, 0, and super 
fragmentation and coherent fragmentation [3| . 

For the spinor Bose-Einstein condensates of bosons 
with spin one, the pairwise interaction Ua is spin- 
dependent. It can be written as IsL Igl ITol lll| 



C/y- = 5{ri-rj){graPm + gqPq) 



(1) 



where Pm and Pq are the projectors to two-body spin- 
states with two-body spin equal to and 2, respectively, 
and gm and gq are the associated strengths. F; is the 
spin operator of particle i, cq = (gm + '^gq)/'^ and C2 = 
{gq — gm)/'i. The term with cq is a central interaction, 
while the term with C2 is spin-dependent. 

When gq gm (e.g., for ®^Rb), we have C2 << co and 
therefore the spin-dependent term is unimportant. In 
this case the single-mode approximation (SMA) works 
very well a common mode 

function (pijc) is obtained by solving the Gross-Pitaevskii 
(CP) equation with the central interaction alone, while 
the small spin-dependent interaction is treated as a per- 
turbation. 

When C2 is not very small, the perturbation treatment 
is poor. In this case, a three-mode approximation (TMA) 
is proposed. Instead of (/)(r), three mode functions 0^(r) 
are introduced, where fi = 1,0, and —1 are associated 
with the three components of spin. These functions obey 



a set of coupled GP eq uations with all the interactions 
taking into account [iJl • 

In the above theories the total spin S of the conden- 
sate is not explicitly considered but only its Z-component 
Sz- However, the interaction (eq.(Q) assures that S is 
exactly a good quantum number. In fact, the eigenstates 
should be S'-dependent but S'^-independent (if external 
magnetic field is not applied) . Therefore, a more per- 
fect theory should preserve S . It is the aim of this 
paper to introduce a ^-conserved description for spinor 
condensates. For this purpose, we need the knowledge 
from few-boson systems. This is discussed in the next 
section. 

2, Few-boson model systems 

Let us first inspect some few-boson 2-dimensional 
model systems with spin one and mass to, and they are 
confined by a harmonic trap with frequency to. When 

hw and \ are used as units of energy and length, re- 
spectively, the Hamiltonian of the TV-boson system reads 



H 



(2) 



Using the method as given in the ref.0,^3 , the Hamil- 
tonian for iV = 2, 3, and 4 were diagonalized. Of course, 
both the orbital angular momentum L and S are good 
quantum numbers. However only the L — states are 
considered here. For a series of eigenstates with a given 
S, let the lowest one be called a first-state. Let the 
energy and wave function of a first-state be denoted as 
Es and '^s ■ When gq + gm = 0.4 is assumed, the 
evolution of Es in accord with f3 = gq/ gm is given in 
Fig. 1. In this figure, there are three groups of curves 
for the three cases = 2, 3, and 4, respectively. In each 
group the curves have distinct S but the same Sz ~ 0. 
When (3=1, the interaction is spin-independent and the 
ground state is degenerate against S. When fi deviates 
from 1, the interaction depends on the spin and, accord- 
ingly, the ground state splits into a band, the ground 
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band. For an example, when N = 4, the ground band 
contains three members with S — 0, 2, and 4. When 
/3 < 1, a pair of boson would be less repulsive if their 
spins are coupled to 2 instead of (because gq < gm). 
Since in ^4 any pair of bosons must couple to 2, this 
state is lower than '52 and 'I'o. On the contrary, when 
P > 1, '^4 is the highest state of the band. Since the 
split of the ground state arises from the character of the 
interaction, it would remain when N is large. In this 
case the spins will have many ways to couple with each 
other to form total spin-states with S ranging from (or 
1) to TV. It implies that , when N is large, the ground 
band would contain many members distinct in their to- 
tal spin-states. To clarify the structure of the ground 
band is another aim of this paper. Incidentally, since 
the members of the band are already orthogonal due to 
their spin structures, their spatial wave functions may be 
alike. 

Let us analyze the wave functions to see the structure 
of the ground band of the few-boson systems. In the 
^-conserved form, the wave function is written as 

^-5 = 51 fs^^^SM (3) 

Xi 

where 6s, m is the total spin-state with the good quantum 
number S, and it is also the i — th basis function of the 
A representation of the permutation group (e.g., refer to 
ref.[15]). fsxi is a function of spatial coordinates be- 
longing to the same representation. The summation over 
i is necessary to assure 'I's to be totally symmetric. The 
summation over A responds to the possibility that some 
total spin-states distinct in A may have the same S (e.g., 
for = 3 and 5 = 1, there are three total spin-states 
^i,{3}i, ^'i,{2i}i, and 6'i,{2i}2 Y^- From eq.®, for a *s, 
we can define the weight of a representation as 

Wx=Y,< fsM\fsx^> (4) 

i 

It turns out, from our numerical results, that all the 
states of the ground band are dominated by A = {N} 
component. E.g., when = 4, the second lowest state is 
\l/2 (cf. Fig.l). For S = 2, there are six total spin-states 
belonging to the {4}, {31}, and {22} representations, re- 
spectively. When P — 0.6 , and g^ + gq varies in the 
range from to 3, our calculation shows that the W|4} of 
^'2 ranges from 1 to 0.996, it implies that the symmetry 
{4} is dominant. The increase of gm + gq causes only a 
very small decrease of W^{4} from the unity. The sym- 
metries other than {4} are found to be important only in 
excited states but not in the ground band. This holds 
also when (3 is given at other values. The underlying rea- 
son is that the A of the spatial wave function must be the 
same as the total spin-state. If A all the bosons 

can not condensate to the same spatial state, at least a 
boson should be spatially excited. This excitation will 
cause an explicit increase of energy. Since the ground 
band is an aggregation of the first-states, each of them 



will do their best to lower the energy. Therefore the 
{N}-symmetry is pursued. Evidently, the above expla- 
nation remains valid not matter how large N is. Thus, 
it is expected that the ground band of a large- TV system 
remains to be dominated by A = {-^}- 

3, Symmetrized total spin-states 

Let us now consider the ground band of the systems 
with many bosons . Based on the above discussion, only 
the case A = {A^} is necessary to be considered. As a first 
step, we have to clarify the number of symmetrized total 
spin-states. For a iV — boson system, let denotes 
a normalized totally symmetric spin-state having only 
Sz conserved, where the superscript j is introduced to 
number this kind of states. For a given Sz, let the total 
number be denoted as Mz{Sz), we have 

MziSz) ^ -{N - Sz ^-^ ) + l (5) 

On the other hand, let 6g denotes those having both 
S and Sz conserved. For a given S, let the total number 
of dg be denoted as M{S), which is related to Mz{Sz) 
as 

M{S) = Mz{S) - Mz{S + 1) = (1 + (-l)^-^)/2 (6) 

From ® we know that M{S) is one if — 5 is even, 
or is zero if N — S is odd. Therefore, the superscript k 
is redundant, and a simpler notation ^5 is used here- 
after. Evidently, ^5 is a simpler notation of 0s.{n}i 
defined before. For the calculation of matrix elements, 
the following formula is useful. 

<es\F,.F,\9s> = < es\J2F^■F,\9s > 

^ SiS+l)-2N 

N (N-l) ^ ' 

Where the first equality holds because 6*5 is totally sym- 
metric, the second equality holds because 

2^F,.r,= (^F,f -^F^ (8) 

i<j i i 

When 6s are used for the spin-states, the spatial wave 
functions must be totally symmetric just as those for the 
scalar bosons. Therefore, in each member of the ground 
band, condensation would occur just as in the scalar sys- 
tem disregarding what S is. Thus, as a reasonable ap- 
proximation, ^'5 that has been given in eq.® is now 
simplified as 

*5-n,iIi0s(r,)0s (9) 

Where the subscript in (f>s responds to the fact that the 
spatial states may depend on S. Incidentally, the func- 
tions 05 (r) with distinct S are in general not mutually 
orthogonal. 
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4, The generalized Gross-Pitaevskii equation 

Let us insert ^ into the three-dimensional 
Schrodinger equation 



(10) 



where H is given in Using the standard variational 
procedure, a generalized Gross-Pitaevskii (GP) equation 
can be derived as 



'2 
where 



(-V^ + r2) + {N- l)gs Ps]^s{r) = Ss (j^sir) (11) 



9s 



co+ < es\Fi.F2\9s > C2 

S{S+l)-2N 



Co 



N {N -1) 



-C2 



(12) 



0s (r) is normalized and isotropic, ps — (t> s {''')* (f' si''') 
and £5 is related to the eigenenergy Es as 

N (N ~1) 

Es^Nes ^— '-gs p's (13) 

where 

P's= J dr(j}s{ryps(l}s{r) (14) 

The generalized GP equation is the same as the usual 
one with only one exception, namely the strength gs is 
not a constant but S— and A'^— dependent. Consequently, 
the bosons in each state of the ground band condensate 
into a specific single-particle states <ps depending on S. 
From ((HJ we know that the band contains totally A^/2+1 
states if N is even with their 5* ranging from 0, 2, 4, 
to N, or contains [N + l)/2 states if N is odd with S 
ranging from 1, 3, 5, to A^. When 13 = 1, the generalized 
GP equation will automatically reduce to the usual GP 
equation. When f3 > 1, C2 is positive, a larger S will lead 
to a larger gs , and therefore a higher level. Thus, the 
upmost level of the ground band has S = N , while the 
ground state has 5 = 0. On the contrary, when /3 < 1, 
the ground state has S = N. This finding coincides with 
previous findings 0|. 

Let us go to realistic cases (three-dimensional) with a 
large N to inspect the structure of the ground band. For 
^^Na atoms, we have g^, — 6.351 x lO^'^v^ a-nd /3 = 1.10 
. For ®'^Rb atoms, we have gm = 2.52 x 10~^y/uj and 
f3 = 0.986. By solving lfTT)l . the ground state energies 
Eg per particle with respect to Nuj^^^ are plotted in Fig. 

2 . In order to see the broadening of the ground band, 
(Es - Es^^J/N with respect to S'/A^ is plotted in Fig. 

3 , where uj = lOOHz is given, the subscript Smin — (if 
N even), or = 1 (if N odd). 

Now, let us evaluate the level-density. When S/N << 
1, eq. (Illll can be solved perturbatively (the perturbation 



is caused by gs — gs^i„), we obtain 

r^r N - S{S + 1) - 

£s-£s„,„ = (^-l)(5s-gs„,J Ps^,„^ C2 —-^ — - Ps^,„ 

(15) 

Thereby the difference Es — Es^^^ can be obtained, 
and we have the level-density 

D{Es) = I — | = l/(|c2|Ps_ (25+1)) 

= S{S + l)/(2\Es-Es„^J{2S+l)) (16) 

where dn is the number of levels in dEs- Since 
\Es — Es^i^l is a very small quantity in the vicinity of 
ESmin I the levels are extremely dense in the vicinity of 
the ground state (if /? > 1) or in the top of the ground 
band (if /3 < 1). This is an interesting point. 

When S/N is not very small, via a numerical evalua- 
tion with a high precision, Es can be represented by a 
parabolic curve 



Es - i?s„„„ ± 



(17) 



where d is a function of N as given in the inset of Fig. 
3 , the positive (negative) sign is for /3 > 1 (< 1). The 
band width Ei,^ = \Em — Es^^^ \ is given in Fig. 4 versus 
Nuj^/^ . The level-density is then 



DiEs) 




N 



Es 



(18) 



They are given in Fig. 5, where uj = lOOHz is assumed. 

If the generalized GP equation is not solved exactly 
but perturbatively (the perturbation is caused by gs — cq 
) , then the resultant energies are identical to those given 
by the SMA^. 

5, Final remarks 

In conclusion, we have derived a generalized GP equa- 
tion based on a 5 -conserved formalism . This equation 
can describe not only the ground state but the whole 
ground band. In the band the states having larger S 
are higher (lower) in Es if /?>(<) 1 . It is recalled 
that S is in fact a good quantum number. The present 
theoretical framework satisfies this requirement, and has 
a very simple formalism to facilitate numerical calcula- 
tions. Furthermore, ^ is the only approximation, which 
is comparable to that used for deriving the well known 
GP equation. Therefore, the reliabilities of these two 
equations are comparable. Thus, this work is an im- 
provement of the SMA and TMA, while at the same time 
avoids perturbative treatment or the difficulty of solving 
coupled equations. 
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FIG. 1: The evolution of Es (in hui ) with /? in few -boson 
systems. Where the lowest two curves, the two in the mid- 
dle, and the upper three are for 2-, 3- and 4-boson systems, 
respectively. S is marked by the Curves. A is introduced to 
shift the curves, A = 0, 1, and 2 for Af = 2, 3, and 4 systems, 
respectively. 




FIG. 2: The ground state energy Eg per particle (in huj) 
versus Noj^^'^. The total spin of the ground state is or 1 
(if ^^Na), or is iV (if **^Rb). In Fig. 2 to 5 the solid (dashed) 
lines are for ^^Na (*^Rb). The inset is just for smaller Nui^^^ 
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FIG. 3: {Es- Es^i„ ) /N (iu liu) versus S/N with a; = WOHz 
. The four sohd (dashed) hues from the lowest to highest 
(highest to lowest) are for iV = 10^ 10^ 10^, and 10^ re- 
spectively. The inset is d versus N. 
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FIG. 5: The level-density D{Es) versus Es - Es^^^ (in tko) 
with u} = lOOHz. The lines from the shortest to longest are 
for TV = 10^ 10^, 10'', and 10^ respectively. 



